Solitons in PT-symmetric nonlinear lattices 
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Existence of localized modes supported by the PT-symmetric nonlinear lattices is reported. The 
system considered reveals unusual properties: unlike other typical dissipative systems it possesses 
families (branches) of solutions, which can be parametrized by the propagation constant; relatively 
narrow localized modes appear to be stable, even when the conservative nonlinear lattice potential 
is absent; finally, the system supports stable multipole solutions. 
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PACS numbers: 42.65.Tg, 42.65.Sf 

Since introduction of the concept of the PT-symmetric 
potentials this subject attracted a great deal of atten- 
tion 0. While the primary interest was devoted to such 
systems in the context of non-Hcrniitian quantum me- 
chanics, recently new applications of the PT-symmetric 
potentials have been found in optics in media with inho- 
mogeneous in space gain and damping, i.e. with properly 
designed imaginary part of the linear refractive index. 
The first experiments reporting the phenomenon are al- 
ready available Q . As soon as the importance of the op- 
tical applications was realized, it became also clear that 
the phenomenon can be studied in the nonlinear context, 
from the point of view of existence of nonlinear local- 
ized modes in linear PT-symmetric potentials It is 
then natural step to address the existence and stability 
of localized modes in nonlinear PT-symmetric potential 
which in optics can be implemented by means of proper 
spatial modulation of nonlinear gain and losses. As an ex- 
ample, such optical systems can be nonlinear waveguides, 
employing concantenated semiconductor optical ampli- 
fier and semiconductor doped two-photonic absorber sec- 
tions (notice that experimental implementation of the lin- 
ear PTsymmetry breaking was reported in Q). 

While it is now known that stable localized [^, 0] and 
moving 0] solitons can exist in conservative purely non- 
linear lattices (see also 0, [l^ for review) the ex- 
istence of stable localized solitons in complex nonlinear 
lattices is still an open problem, since up to now only 
periodic waves were found to be stable in such struc- 
tures 111. The elucidation of stable localized solitons 



in PT- symmetric nonlinear lattices is therefore a central 
goal of this Letter. 

We describe the propagation of laser radiation along 
the ^-axis of the medium with periodic transverse mod- 
ulation of cubic nonlinearity and nonlinear gain with the 
complex nonlinear Schrodinger (NLS) equation for the 
dimensionless light field amplitude q: 

m = Alvv - kP? - [V{V) + ^W{rj)] \q\'q (1) 



gitudinal coordinates, respectively. The functions V{ri) 
and W{r]) describe transverse periodic modulations of the 
conservative and dissipative parts of the nonlinearity and 
are assumed to satisfy the T^Tsymmetry relations. We 
further assume that conservative and dissipative parts 
of nonlinearity have the same period tt, i.e. V{r]) = 
V{-ri) = V{r] + tt) and W{r)) = -W{-ri) = W{7] + tt). 
These functions will be considered bounded with ar and 
CTi being the maxima of V and W, respectively. 

We are interested in stationary localized solu- 
tions, which can be searched in the form q{T], ^) = 
u{ri)e^^'^^^^^'^^ = [wrijfj + iwi{r])]e^''^ where u and 9 are 
the amplitude and phase of the mode. Eq. ^ can be 
rewritten in the hydrodynamic form 
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-bu+[l + V{v)y^^ 



0, jr, = ~2W{ti)u^ (2) 



where we have introduced the " current density" given by 
j = u^dO/dr]. 

Starting with general properties of stationary localized 
solutions, we notice that it follows from ^ that such 
solutions can exist only for 6 > and their asymptot- 
ical behavior at ±oo is given by w ^ g±V2brj 
j ~ g±4%/2hi7^ J ^j^g current density is localized much 
stronger than the field. One easily finds that the field u 
can become zero only in the points where the current den- 
sity j is zero, as well. Furthermore, closely following the 
approach described in [§], one obtains that b = 0{U^) 



where U = Ur + Ui, with Ur.i = 



iiri)dr], is the 



where rj and ^ arc the normalized transverse and Ion- 



total energy flow of the beam. This implies that in the 
limit of small intensity {U — )■ 0) we have U ^ \fb, and 
respectively Umaa; — > 0. Following one can obtain 
the relation h < u'^^^{l + ai). Thus the existence of the 
solutions with Umax implies 6 — >■ 0. 

The energy balance in stationary solutions follows from 
^ and reads W{r])u'^{r])dr] = 0. Since we consider 
odd functions W{r]), this condition can be satisfied by 
any even function u{r]). In other words, unlike this hap- 
pens in dissipative systems of a general kind the 
requirement of balance between losses and gain in our 
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system does not introduce a constraint selecting only one 
possible mode (i.e. the propagation constant is not deter- 
mined by the balance between losses and gain). Thus, in 
terms of the existence of branches of solutions the proper- 
ties of Eq. ([T]) resemble the properties of the conservative 
NLS where the propagation constant b is determined by 
the energy flow U and a continuous family of solutions 
exist. This fact is illustrated in Fig. [T] Hereafter in all 
numerical simulations we use 

V{r]) ^ arC0s{2ri) and VK(r/) = -cr, sin(2?/), (3) 

where and ai are the modulation depths of the conser- 
vative and dissipative lattices. In Fig. [1] (a) we observe 
that increase of b results in monotonic growth of the soli- 
ton energy flow U and the contraction of light in a single 
channel of nonlinear lattice. Now, however the energy 
of the soliton is distributed between real and imaginary 
components of the field, as it is shown in Fig. [T](b). The 
ratio Ui/Ur of energy fiows concentrated in imaginary 
and real parts of the field takes on maximal value at in- 
termediate b values and diminishes at 6 — )- and 6 — >■ cx). 
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FIG. 1: (Color online) (a) Energy flow versus propagation 
constant for fundamental solitons at (T; = (curve 1), 1.4 
(curve 2), 1.8 (curve 3) for Or = 1. (b) Ui/Ur versus propa- 
gation constant at ai — 1.8, ar — 1. 



Let us now turn to more detailed study of the men- 
tioned limits of the propagation constant. First of all. 
Fig. [T] (a) supports the above estimate U ^ b^^^, also 
illustrating that in the limit 6—5-0 the energy fiow very 
weakly depends on the amplitude of dissipative part of 
potential (the three lines arc indistinguishable on the 
scale of the picture). Indeed, in this case the lattice 
period becomes small in comparison with the width of 
smoothly modulated soliton and one can perform the 
averaging procedure 0]. For the model ^ the solu- 
tion of ([ij can be found in the form <z(?7,'^) ~ Qiv^O + 
A{7], ^) cos(2?y) + iB{7], ^) sin(2r/), where Q, A, and B are 
the functions slowly varying on the scale tt. Substitu- 
tion of this ansatz in Eq. ([1]) yields A = (Tr\Q\'^Q/2 and 
B = — (TilQpQ, and cubic-quintic NLS equation for the 
field Q: 



where x = ^(^c^ — erf) This equation does not contain 
any imaginary part - the consequence of the opposite 
parities of real and imaginary components of the nonlin- 
earity modulations. The solitonic solution of (j4]) which 



exists at &x > — 1/8 is well known (see e.g. jl2l |) 



Q = 2Vbe'''^ \l + VTT8x6cosh(2\/2677) 



-1/2 



(5) 
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iQi + 7^Qnn + \Q\Q+oX\QrQ = 0- 



(4) 



This solution is reduced to the conventional NLS soliton 
in the limit 5 — )■ 0, revealing weak dependence of the 
soliton on the parameter x, what explains convergence 
of all branches in Fig. [T] (a) at 6 0. 

The profiles of the simplest fundamental soliton solu- 
tions of Eq. ([1]), i.e. the solitons belonging to the low- 
est branch (see also Fig. |3] (a), below), are shown in 
Fig. [2] (the phase of the solution is fixed by the condition 
6{0) = 0, it however can be changed due to the phase 
invariance of the complex NLS equation). The centers of 
such solitons reside in the point where conservative part 
of nonlinearity takes on the maximal value, while dissi- 
pative part of nonlinearity is zero. Due to the fact that 
left wing of soliton resides in the domain with nonlin- 
ear losses, while its right wing is subjected to nonlinear 
gain the solitons are characterized by the anti-symmetric 
imaginary parts of the field [Figs. [JJa) and (b)] indicat- 
ing on tilted phase fronts and the existence of internal 
currents directed into the domain with losses [Fig.[2](d)]. 

Returning to the simple approximation ([5]) we also ob- 
serve that at fixed b and ar it suggests the existence of 
the upper limit, < cr,"^^, of the strength of the dis- 
sipative term ai, for which localized dissipative solitons 
exist. This is indeed confirmed numerically in Fig. 
[notice that the simple estimate for this upper limit 
^(upp) ^ 2^2 _|_ 3/4 gives for the parameters of Fig. [3Ja) 
a'^PP Ri 1.22 while the numerical value is a"^''' w 1.62]. 
The growth of ai results in the monotonic increase of 
the imaginary part of the field [c.f. also Figs. [2][a) and 
(b)] accompanied by a considerable increase of current 
density [Fig. HJ^d)]. The energy flow increases with ai 
[Fig. ISlJa), red curve] until the tangential line to U{ai) 
becomes vertical. Apparently, there exists another up- 
per branch of solutions joining with the lower branch in 
the point ai = ct"^^ [Fig. Elja), black curve] for which 
the energy flow is a monotonically decreasing function of 
ai. The solitons belonging to this branch are character- 
ized by a double-hump field modulus profile [Fig. [2jc)]. 
When ai decreases the real part of the solutions decays 
and only imaginary survives. The later is asymmetric 
and its maximum and minimum are located in a single 
period of V{ri) [this tendency is visible in Fig.[ljc)]. The 
solitons from upper branch in Fig. [3] (a) are unstable. 
Besides these simplest branches one can find a variety 
of soliton families with more complicated internal phase 
distributions, but we do not discuss them here because 
they are usually unstable. 

One of the most important results of this Letter is that 
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FIG. 2: (Color online) The profiles of fundamental solitons 
from lower branch at ai = 0.5 (a) and 1.58 (b). (c) The profile 
of soliton from higher branch at ai = 0.3. (d) The current 
density for the fundamental solitons from lower branch at 
CTi — 0.5 (curve 1) and 1.58 (curve 2). In all cases b = 1 and 
CTr = 0.5 



fundamental solitons can be stable despite the fact that 
the system ([1]) is characterized by the presence of do- 
mains where only losses or gain are acting. The outcome 
of stability analysis is presented in Figs. [3] (b)-(d). The 
fundamental solitons are stable for below certain crit- 
ical value erf [see Fig. [3] (d) for a typical dependence 
of the perturbation growth rate on ct^]. Notice that for 
ai > af the growth rate 5r increases until one reaches 
the border of existence domain Ui — a^^^. For fixed <7r 
the stability domain on the plane {b, ai) is rather complex 
[Fig. 131(c)]. At ar = 0.5 both af and cr"^^ increase as 
— > indicating on soliton stability in a very broad range 
of amplitudes of gain modulation. For sufficiently large 
b values the domains of existence and stability monoton- 
ically expand with b. A similar situation is encountered 
for other values of ar. The increase of the depth of modu- 
lation of conservative nonlinearity ar at fixed b results in 
considerable expansion of existence domain on the plane 
{ar, ai). 

However, especially interesting situation occurs at 
ar = 0. In this case, there is no modulation of conser- 
vative nonlinearity at all, but our analysis still predicts 
stability of fundamental solitons between two red lines 
in Fig. ^c) (for b > b'^^ ~ 1.05 the solitons are stable 
for < ai < af). This fact is really remarkable tak- 
ing into account that now the symmetric conservative 
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FIG. 3; (Color online) (a) The energy flow vs ai for lower 
(red curve) and upper (black curve) branches of fundamental 
solitons at b = 1, ar ~ 0.5. Circles correspond to solitons 
shown in Figs. [Ha)-(c). Domains of existence and stability 
on the plane {b,ai) for fundamental solitons at ar — 0.5 (b) 
and ar — (c). (d) The perturbation growth rate vs ai at 
ar = 0.5. 

nonlinearity providing the restoring force in the case of 
slight displacements of soliton center from the equilib- 
rium position 77 = 0, is absent. We observe that the loss 
of stability occurs at the soliton width, which is compa- 
rable with the characteristic scale of the lattice, i.e. to 
the half-period 7r/2. Thus, the modulation of conserva- 
tive nonlinearity is not a necessary ingredient for soliton 
stability, although it can change considerably stability 
properties of low-power solitons with 6 — )■ 0. 

In addition to the fundamental solitons we found a 
variety of multi-hump states whose humps reside on dif- 
ferent maxima of conservative nonlinear lattice V . The 
representative examples of such states that in the limit 
CTi — >■ transform into conventional dipole and tripole 
solitons are shown in Fig. |4l In such solitons the real 
part of the field (dominating at ai — > 0) changes its sign 
between neighboring maxima of V . The current density 
in such states is characterized by n (n is the number of 
peaks in field modulus) negative spikes in the vicinity 
of maxima of V. Analogs of solitons with in-phase field 
peaks were obtained too, but they all are unstable. Like 
fundamental solitons, multipole states are parameterized 
by the propagation constant b. For a given ar there exist 
a cutoff on b below which multipole solitons do not exist, 
while increase of b results in growth of energy flow. 

Increase of gain-loss modulation ai also causes increase 
of U and fraction of power concentrated in imaginary part 
of the of multipole soliton [c.f. Figs. IDJa) and (b)], but 
such solitons can be found only at ai < af^^ [in Fig. [DJa) 
we show only the lower branch of dipole solitons although 
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FIG. 4: (Color online) Dipole solitons at (a) ai = 0.5, (b) 
CTi = 2.67, and tripole solitons at (c) ai = 1. In all cases 
ar = 1, 6 = 2 

upper unstable branch can be found too]. Linear stabil- 
ity analysis predicts stability of the niultipole solitons at 
CTi < erf as shown in Fig. El^c). This domain of stability 
gradually broadens with increase of the depth of mod- 
ulation of conservative nonlinearity ar [Fig. [5] (b)]. In 
contrast to fundamental solitons multipole solitons can 
be stable only if propagation constant is sufficiently large. 
This critical value of propagation constant increases with 
decrease of Cr- This is because multipole solitons may 
exist only if conservative nonlinearity is modulated and 
when this modulation is sufficient for compensation of 
repulsive forces acting between neighboring poles. In- 
crease of the number of poles in solitons does not result 
in dramatic modifications of existence domain [0,(t"^^] 
but domain of stability [0, erf] shrinks with n. 

To conclude, we have reported a set of stable local- 
ized solutions supported by T'T-symnietric nonlinear lat- 
tices. The system considered reveals a number of unusual 
properties. First, although it is dissipative and the bal- 
ance between the gain and losses must be satisfied, it 
possesses families (branches) of solutions, which can be 
parametrized by the propagation constant 6, in contrast 
to other typical dissipative systems. Second, the modes, 
whose width is smaller than the lattice half-period appear 
to be remarkably stable, even when the conservative non- 
linear potential is absent. Finally, the system supports 
stable multipole solutions. 

The work of FKA and VVK was supported by the 
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FIG. 5: (Color online) (a) Energy flow vs nonlinear gain for 
dipole solitons at a,- — 1, b — 2. Only lower branch is shown. 
Circles correspond to profiles shown in Figs.[4][a) and (b). (b) 
Domains of existence and stability on the plane {ar,(Ji) for 
dipole solitons at 6 = 2. (c) Real part of perturbation growth 
rate vs the gain parameter at ar- 
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